Towards the electron EDM search. Theoretical study of PbF 
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We report ab initio relativistic correlation calculations of potential curves and spectroscopic con- 
stants for four lowest-lying electronic states of the lead monofluoride. We also calculated parameters 
of the spin-rotational Hamiltonian for the ground 2 n 1/ / 2 and the first excited A 2 E]^ 2 states including 
P. T-odd and P-odd terms. In particular, we have obtained hyperfine constants of the 207 Pb nucleus. 



For the II1/2 state A 



-6859.6 MHz, A« 



9726.9 MHz and for the A 2 E+ /2 A A 



1720.8 MHz, 



Am = 3073.3 MHz. Our values of the ground state hyperfine constants are in good agreement with 
the previous theoretical studies. We discuss and explain seeming disagreement in the sign of the con- 
stant A± with the recent experimental data. The effective electric field on the electron E e g, which 
is important for the planned experiment to search for the electric dipole moment of the electron, is 
found to be 3.3 x 10 10 V/cm. 
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The interest in theoretical study of 207 PbF is caused 
by the planned experiment to search for the simultaneous 
violation of the time-reversal invariance (T) and space 
parity (P) The latter may be caused either by the 

P, T-odd electron-nucleus interaction or by the electric 
dipole moment of the electron (eEDM, or d c ), which can 
exist also only due to fundamental P, T-odd interactions 
[f|. The possible experiments to search for P, T-odd ef- 
fects in atoms and molecules were discussed in several 

books and reviews 0-112. 

Experiment to search for cEDM with PbF molecule 
is currently under preparation by the group of Shafcr- 
Ray Q . PbF molecule has one of the strongest effective 
electric field E e g on the valence electron, which deter- 
mines experimentally measured frequency shifts [l], 
This experiment can test "new physics" beyond the Stan- 
dard Model including different supersymmetric models. 
Many such models predict eEDM in the range between 
10 -28 — 10~ 29 e ■ cm (e is the electric charge of the elec- 
tron) , within the reach of the new generation of experi- 
ments mm- 

Recently McRaven et al. 0] have measured hyperfine 
structure parameters for the ground and the first excited 
states of 207 PbF and found them to be in disagreement 
with previous predictions 0, Q ■ Taking into account that 
theoretical studies of the lead monofluoride took place 
decades ago within rather simple theoretical models we 
recalculate spectroscopic constants and parameters of the 
spin-rotational Hamiltonian, H sr for PbF molecule. In 
particular, we focus on the calculations of the effective 
electric field on the electron E e f{ and the hyperfine struc- 
ture constants for the lead nucleus. Our new results are 
in reasonable agreement with previous calculations 0, Q ■ 
We discuss the discrepancy between the theory and ex- 
periment and argue that it is caused by the incorrect 
phase factor in the spin-rotational wave function, rather 
than in the difference in the sign of the hyperfine constant 
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The Hamiltonian H sr can be written as [2 



H sr = B'J 2 + AS' ■ J + S'AI + (1) 
+^ S'GB — Dn ■ E + (2) 

+W P K P n x S' • I + (W P , T Kp,T + W d d e )S' ■ n, (3) 

where the first line corresponds to the rotational struc- 
ture with Si-doubling and the hyperfine interaction of the 
effective electron spin S',5" = \, with the spin of the 
nucleus I; H is the rotational constant, J is the electron- 
rotational angular momentum, A is the ^-doubling con- 
stant, A is the hyperfine structure tensor. The second 
line describes interaction of the molecule with the exter- 
nal fields B and E, n is the unit vector along the molecu- 
lar axis, C, directed from Pb to F; G is the G-factor ten- 
sor. The last line corresponds to the interaction with the 
T-forbidden anapole moment of the nucleus, the P, T-odd 
weak interactions of the electrons with the heavy nucleus 
and the interaction of eEDM with the internal electric 
field of the molecule; up and kp,t are the anapole con- 
stant and the P, T-odd neutral current constant for Pb 
nucleus. 

Parameters W% in ([3]) can not be measured experimen- 
tally and have to be found from the molecular electronic 
structure calculations. Their accurate knowledge is cru- 
cial for extracting constants up, np,r, and d e from the 
experiment. Expressions for calculating parameters Wi 
can be found in Refs.[l|, 0, [ll|. In the literature the ef- 
fective electric field seen by the unpaired electron is usu- 
ally used instead of the parameter Wd- It is defined as: 
E cS = W d \to\, where ft = (*fi|S' • n|*n) = ±1/2 and 
\1/q is the electronic wavefunction for the 2 n!/ 2 or 2 Sy 2 
states. 

We used the generalized relativistic effective core po- 
tential (GRECP) to simulate the interaction of va- 
lence electrons (4 electrons of Pb) with the explicitly ex- 
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eluded Is to 4/ electrons of Pb (68 core electrons). In 
addition the 5 s, 5p and 5 c? shells of Pb (treated explicitly 
with the used GRECP) were frozen. All the electrons of 
fluorine were treated explicitly. The resulting configura- 
tions were: 6s 2 6p 2 for Pb and ls 2 2s 2 2p 5 for F. For Pb 
we used generalized correlation atomic basis set [l5, 16 1 



pJ5sl6pl2o!9/)/[5s7p<W2/] a nd GRECP from the paper 
For F we used the ANO-I (Us9p4d3f)/\4s3p2dlf} 
atomic basis from the MOLCAS 4.1 library [17j . 

By means of the complete active space self-consistent 
field method the molecular orbitals were received. Molec- 
ular symmetry was treated by using the C 2v point group 
classification scheme. The next step was to apply the 
spin-orbit direct configuration interaction (SODCI) ap- 



proach [181-1201] (modified in [2l( to account for spin-orbit 
interaction in configuration selection procedures). The 
main idea of this method is to use selected singly and 
doubly excited configurations with respect to some mul- 
ticonfigurational reference states within the second order 
perturbation theory. After that the SODCI approach is 
employed on the space of all the selected configurations. 

In the present work we calculated four lowest states, 
which belonged to the configurations [...]<t 2 tt 1 / 2 ( 2 ^-i/2)i 
[..V*-3/a ( 2 n 3/2 ), [...]aV (A 2 S+ /2 ) and [...]aV 
(B 2 £^ 2 ). Here a is mainly formed from the 6 s basis 
function of Pb, tt\/ 2 and n 3 / 2 are mainly made from 6p x 
and 6p y functions of Pb. The orbitals a' and a" are 
formed from the 6p z and 7s functions of Pb, respectively. 

One of the main goals of the present paper was to 
draw potential curves for the ground and three lowest 
excited states. Therefore, we made calculations for sev- 
eral internuclear distances (from 3.8 a.u. to 4.2 a.u.) near 
the ground state equilibrium distance (i? c = 2.0575 A « 
4.0 a.u.). 

Equilibrium internuclear distances, vibrational con- 
stants, T c and dipole moments for four lowest states of 
PbF molecule can be found in Table U (the most pre- 
cise data, with 0.00005 treshold, are pointed out for the 
D(R)). The PbF properties could be divided into two 
groups: valence and core ones. The parallel and per- 
pendicular G-factor components, the electronic dipole 
molecule moment D(R) or fi-doubling constant could be 
related to the valence properties. However, such prop- 
erties as E e g, A±, Au, anapole moment constant Wp 
or T-odd interaction constant Wp,r could be referred to 
the core ones. Mainly all these properties were calcu- 
lated for the ground and first excited states at the point 
4.0 a.u., which is close to the equilibrium distance. In 
the core region one should make a four-component rela- 
tivistic calculation. However, we have performed calcula- 
tions on two-component orbitals by means of pseudopo- 
tential and, therefore, have to apply a restore procedure 
for four-component orbitals from two-component ones 
at the next stage. For this purpose the nonvari ational 
one-centre restoration method developed in 



is employed. One can find the calculated parameters 
of the spin- rotational Hamiltonian, H sr , in Table IlLTl for 
the 2 ni/2 and A 2 E^ 2 state, respectively (where (a) and 
(b) of [l| correspond to the minimal and maximal spin- 
orbital mixing) . Finally, in Table |H] the matrix element 
values of Sl-doubling constant A divided by 2/ 2MB 2 
(A/2B'(R)) and electronic dipole moments are listed for 
five investigated points and for two electronic states of 
the molecule. 



TABLE I: Internuclear distances, vibrational constants, T c 
and dipole moments for four lowest states of PbF molecule 







-Re (A) 


w e (cm x ) 




D(R) (D) 


*n 1/2 


m 


2.06 


489 




4.62 




Experiment 


2.0575 


502.7 









This work 


2.08 


493.6 





4.26 




Experiment 


2.0342 


528.7 


8263.5 




This work 


2.07 


522.3 


8258.9 






[2] 


2.18 


418 




5.5 


Experiment 


2.1597 


394.7 


22556.5 






This work 


2.18 


388.5 


23380 


2.51 




Experiment 


1.97 


605.7 


35644.4 






This work 


1.99 


597.8 


36583.8 
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According to the data pointed out in Table U T c values 
are in good agreement with the experimental data. The 
obtained equilibrium distance, R e , have slightly higher 
value than in the experiment, whereas the presently ob- 
tained vibrational constants are in much better agree- 
ment with the experiment comparing to those from 
Note that the dipole moment D(R) value for the ground 
state are nearly the same to those obtained in Q , whereas 
they are decreased approximately two times for the ex- 
cited state. All data in Table [TTI obtained for D(R) and 
^-doubling constant have been calculated with the 0.0005 
threshold. But we should note that we have computed 
these parameters for point 4.0 a.u. more accurately with 
0.00005 threshold. For f2-doubling these values are -0.361 
and 1.46 for the 2 n!/ 2 and A 2 E^ 2 states, respectively. 

Let us thoroughly examine the data obtained in Ta- 
ble IIIII First of all we would like to discuss the results 
for the ground state. The first point to note is that the 
Ay value is in the fault bounds, meantime the G|| value 
is increased in two times according to Q. The effective 
field constant Wd obtained is the best datum to-date, 
the anapole moment constant Wp and T-odd interac- 
tion constant Wp t T increased their values nearly by 50%. 
The most interesting result concerns the perpendicular 
hyperfine structure constant Aj_. The sign of the ob- 
tained value exactly agree with those in previous works 
0, 0] ■ Since the constants of the H sr such as Aj_ , G±_ , 
^-doubling are determined by the off-diagonal matrix el- 
ement, it is important to fix the phase factors for the 
wave functions to get the consistent signs of those ma- 
trix elements. In our calculations the phase factor is the 
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TABLE II: Matrix element values of f2 - doubling constant divided to 2B'(R) and electronic dipole moments (ID = 0.3934 
a.u.) for the ground and the second excited states of PbF molecule 



H SI constants 


<1<|A/2B'|4<> 


State 


a n 1/2 


A ^1/2 


Internuclear distance 


3.8 


3.9 


4.0 


4.1 


4.2 


3.8 


3.9 


4.0 


4.1 


4.2 


[21 






-0.347 
















i 






-0.303 










1.49 






This work 


-0.412 


-0.342 


-0.369 


-0.391 


-0.411 


1.49 


1.42 


1.45 


1.47 


1.49 




D(R) (D) 


State 


a n 1/a 


A ^1/2 


Internuclear distance 


3.8 


3.9 


4.0 


4.1 


4.2 


3.8 


3.9 


4.0 


4.1 


4.2 


[2] 






4.62 










5.5 






This work 


4.27 


4.264 


4.28 


4.266 


4.255 


2.32 


2.42 


2.55 


2.67 


2.79 



same as in eq. (14) of [![. 

In the experimental work [4] of the Shafer-Ray group 
the obtained A± values are nearly the same by absolute 
value but have the opposite sign for both states. How- 
ever, it was noted in [4| that the energy level structure 
is determined by the expression, in which Aj_ is multi- 
plied by p, where p — ±1 is the parity of the electronic- 
rotational state. Therefore, it was suggested that the 
parity rather than sign of A± has been incorrectly as- 
signed in previous works. However, we believe that the 
changes in parity or sign of the perpendicular component 
are not necessary. Some corrections for the energy level 
expression are required instead. As we found, (~1) J + 1 / 2 
should be replaced by (— 1) J_1 / 2 in eqs. (9) and (23) 
of [H . This leads to change of sign in front of pA± in 
eq. (24) of [l[ for the energy level structure. 

Let us consider the symmetry properties of wavefunc- 
tion of a diatomic molecule at the space inversion opera- 
tion in more details. Let x, y, z be a laboratory axes, and 
the £, 77, £ body frame axes in which the molecule is fixed. 
The a and (3 (7 = 0) Euler angles are defined in such 
a way that the £ axis coincides with the molecule axis. 
The vibrational wavefunction is not important here, and 
electronic-rotational wavefunction within the adiabatic 
approximation in the Hund case c is written as 



(• • • , L, Va, Ca, a a ,...,R) 9 J Mj n(a, P) 



(4) 



where Mj is projection of the total angular momen- 
tum, J, on the z axis, subscript a enumerate electrons, 
a a = ±1 is the spin variable which corresponds to 
the projection of the spin of the a-th electron on the 
£ axis, R is the internuclear distance. ^ fl (a, (3) — 

■\J 2J 4 ^ T 1 D J M J q(Q) j3, 7 = 0) is rotational wavefunction. 
In the laboratory system of coordinates the inversion 
operation, P, is 

P^n (■■■, x a , y a , z a ,a a ..., R, a, /3) 6 J Mjn {a, (3) = 
(. . . , -x a , -y a , -z a , a a . . . , R, a + 7T,7r - j3) 
xOfijAa + K^- f3), 



where = C/*^, U — expiifijf) ■ exp(iaJ^) is the 
unitary operator. Then one can obtain that the inversion, 

P = U PU + , in the body frame system of coordinates is 



■N 



*n (• • • ) -£<», Va,(a, -CT , ...,R) 

x0Mj,n(a + TT,TT- f3), 



(5) 



where i = y— T, N is the number of electrons. Note that 
the above transformation of the electronic wavefunction 
(corresponding to the reflection in the 77 — £ plane) gives 
the electronic wavefunction with the opposite sign of f2, 
which has the same adiabatic potential. In the case of 
Q = electronic wavefunction whether not changed or 
is multiplied by the —1 factor. The rotational wavefunc- 
tions are also transformed into each other with changing 
sign of ft |3: 

The total electronic- rotational wavefunction ((4]), except 
the case O = 0, has no definite parity. In the case 
n two wavefunctions (j4| with different signs of fi are 
transformed into each other, and we can form from them 
positive (p — 1) and negative (p = —1) wavefunctions. 
Let us construct the electronic-rotational wavefunction 
with definite parity for the case £1 = ±1/2 of our in- 
terest, and, in particular, for the 2 n 1 / 2 and 2 E + states. 
Since the closed electronic shell is invariant under reflec- 
tion, consider transformation only of the valence spinor 
<^(£aj Va, Ca> &)> where ui = ±1/2 is projection of the to- 
tal angular momentum, j, of the valence electron on the 
C axis. One can verify that for phase factors accepted by 
eq. (14) of i 

^l/2(-?a,?7a,Ca,-<7) = i(f e _ 1 j% (£ a , X] a , ( a , d) (7) 

and hence, accounting for ([6]) and 0, the electronic- 
rotational wavefunction | J, Mj,p) with parity p becomes 



\J,Mj,p) 



*l/2^M J ,l/2 + (' 



vJ-1/2 



1/2" Mj, -1/2 



(8) 
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Eq. (23) of [1] is reduced to © by replacing (-l) J+1 /2 

on i-iy- 1 / 2 . 

The parallel hypcrfine structure components is in 
agreement with |J|. The other data for the first excited 
A 2 E^ 2 state, listed in the Table Hill were obtained for 
the first time. 

Finally, we would like to note that we have performed 
the most accurate to date electronic structure calculation 
of the lowest states of PbF and, as a result, have obtained 
the most reliable data for the spin-rotational Hamilto- 
nian constants. In this paper we have proposed solution 
of the problem of the sign mismatching in the perpen- 
dicular hyperfine structure constant. The new computed 
data for the spectroscopic constants better conform to 
the experimental results at present. The calculated ef- 
fective electric field on the electron is increased in the 
present study nearly by 10% compared to that in 0j. 

K.B. and A. P. are grateful to the Ministry of Education 
and Science of Russian Federation (Program for Develop- 
ment of Scientific Potential of High School) for the Grant 
No. 2.1.1/1136. A. T. and A.P. were supported by the 
RFBR grant 09-03-01034 



[1]M. G. Kozlov, V. I. Fomichev, Yu. Yu. Dmitriev, L. N. 

Labzovsky, and A. V. Titov, J. Phys. B 20, 4939 (1987). 
[2] . Yu. Yu. Dmitriev, Yu. G. Khait, M. G. Kozlov, L. N. 

Labzovsky, A. O. Mitrushenkov, A. V. Shtoff, and A. V. 

Titov, Phys. Lett. A 167, 280 (1992). 
[3]N. E. Shafer-Ray, Phys. Rev. A 73, 034102 (2006). 
[4]C. P. McRaven, P. Sivakumar, and N. E. Shafer-Ray, 

Phys. Rev. A 78, 054502(R) (2008), erratum: Phys. Rev. 

A 80, 029902(E) (2009). 
[5]L. D. Landau, Sov. Phys.-JETP 5, 336 (1957). 
[6] M. Kozlov and L. Labzowsky, J. Phys. B 28, 1933 (1995). 
[7] I. B. Khriplovich and S. K. Lamoreaux, CP Violation 

without Strangeness. The Electric Dipole Moments of 

Particles, Atoms, and Molecules (Springer- Verlag, Berlin, 

1997). 

[8]E. D. Commins, Adv. At. Mol. Opt. Phys. 40, 1 (1998). 
[9] J. S. M. Ginges and V. V. Flambaum, Phys. Rep. 397, 



63 (2004). 

[10] J. Erler and M. J. Ramsey-Musolf, Prog. Part. Nucl. Phys 
54, 351 (2005). 

[11] A. V. Titov, N. S. Mosyagin, A. N. Petrov, T. A. Isaev, 
and D. P. DeMille, Progr. Theor. Chem. Phys. B 15, 253 
(2006). 

[12] M. Raidal et al., Eur. Phys. J. C 57, 13 (2008), arXiv: 
0801.1826. 

[13] E. N. Fortson, P. Sandars, and S. Barr, Physics Today 

56(6), 33 (2003). 
[14] A. V. Titov and N. S. Mosyagin, Int. J. Quantum Chem. 

71, 359 (1999). 

[15] N. S. Mosyagin, E. Eliav, A. V. Titov, and U. Kaldor, J. 
Phys. B 33, 667 (2000). 

[16] T. A. Isaev, N. S. Mosyagin, M. G. Kozlov, A. V. Titov, 
E. Eliav, and U. Kaldor, J. Phys. B 33, 5139 (2000). 

[17] K. Andersson, M. R. A. Blomberg, M. P. Fiilscher, 
G. Karlstrom, R. Lindh, P.-A. Malmqvist, P. Neogrady, 
J. Olsen, B. O. Roos, A. J. Sadlej, et al. (1999), quantum- 
chemical program package "molcas", Version 4.1. 

[18] R. J. Buenker and S. D. Peyerimhoff, Theor. Chim. Acta 
35, 33 (1974). 

[19] R. J. Buenker and S. Krebs, in Recent Advances in Mul- 

tireference Methods, edited by K. Hirao (World Scientific, 

Singapore, 1999), pp. 1-29. 
[20] A. B. Alekseyev, H.-P. Liebermann, and R. J. Buenker, in 

Recent Advances in Relativistic Molecular Theory, edited 

by K. Hirao and Y. Ishikawa (World Scientific, Singapore, 

2004), pp. 65-105. 
[21] A. V. Titov, N. S. Mosyagin, A. B. Alekseyev, and R. J. 

Buenker, Int. J. Quantum Chem. 81, 409 (2001). 
[22] A. V. Titov, Ph.D. thesis, St .-Petersburg (Leningrad) 

State University, Russia (1985). 
[23] A. V. Titov, Int. J. Quantum Chem. 57, 453 (1996). 
[24] A. V. Titov, N. S. Mosyagin, and V. F. Ezhov, Phys. Rev. 

Lett. 77, 5346 (1996). 
[25] A. N. Petrov, N. S. Mosyagin, T. A. Isaev, A. V. Titov, 

V. F. Ezhov, E. Eliav, and U. Kaldor, Phys. Rev. Lett. 

88, 073001 (2002). 
[26] A. N. Petrov, A. V. Titov, T. A. Isaev, N. S. Mosyagin, 

and D. P. DeMille, Phys. Rev. A 72, 022505 (2005). 
[27] L. D. Landau and E. M. Lifshitz, Quantum mechanics 

(Pergamon, Oxford, 1977), 3rd ed. 
[28] E. R. Meyer and J. L. Bohn, Phys. 

Rev. A 78, 010502(R) (2008), URL 

http : //link . aps . org/abstract/PRA/v78/e010502 



TABLE III: Parameters of the spin-rotational Hamiltonian _ff sr for 2 iI 1 / 2 and A 2 E+, 2 states of PbF molecule 



State 


Reference 


A ± (MHz) 


A l} (MHz) 




G|| 


W d (10 al, )Hz/e cm 


Wp.t (KHz) 


Wp(KHz) 


2 n 1/2 


m 
i 


-8990 
7200 ± 150 


10990 
10300 ± 800 


-0.326 


0.040 


1.4 

1.5 
1.0 
1.8 
1.6 


55 


-0.72 




[28] 
[1] (a) 

[1] (b) 
This work 


-7460 
-8240 
-6859.6 


8690 
9550 
9726.9 


-0.269 
-0.438 
-0.319 


0.034 
0.114 
0.082 


51 
99 
74.95 


-0.65 
-1.25 
-0.99 


A ^1/2 


i 

This work 


-1200 ±300 
1720.8 


3000 ± 2500 
3073.3 


2.417 


1.92 


2.5 


123.7 


-1.588 



